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PRACTICAL CONSIDERATIONS FOR ABEL INVERTING OF 
PHOTOGRAPHIC DATA WITH APPLICATION TO THE ANALYSIS 
OF A 15-kW WALL-STABILIZED ARC-LIGHT SOURCE 

By Walter L. Snow 
Langley Research Center 

SUMMARY 

A wall -stabilized arc source is used to study the practical aspects of photographic 
Abel inversion. The convergence -free problem is set up in terms of energy variables and 
in the limit, the Abel transform derived as it applies to the spectroscopic problem. The 
detailed inversion of the transform makes explicit the implied boundary conditions in the 
solution and also displays two alternative forms to which numerical analysis is usually 
applied. 

The stop is identified which determines the spatial resolution in both photoelectric 
and photographic modes. An expression for the effective spatial resolution is derived in 
terms of the aperture and field stops of the system with due attention to depth of field 
and optical resolution. An upper bound for the relative volume error incurred because 
of converging optics is obtained in terms of geometrical dimensions of the system. 

Temperature profiles were derived from absolute intensity measurements of con- 
tinuum and atomic line radiation. Measurements relate to atmospheric pressure argon 
with maximum flow rate of 3 x 10"5 m3/sec (4 ft^/hr) at root-mean-square power levels 
between 2 and 15 kW. Equilibrium plasma at 8000 K to 12 000 K was generated under 
these conditions. Theoretical and experimental estimates of optical depth for 4300 Arl 
are contrasted. Stark width and shift data for 4300 Arl were used as self-consistency 
checks for the axial temperatures derived from ’’side on" measurements. 

INTRODUCTION 

The popularity of Abel inversion for analysis of cylindrically symmetric phenomena 
is evidenced by the number of references to it in recent spectroscopic literature. Unfor- 
tunately, for one reason or another, the practical aspects can only be gleaned with diffi- 
culty by anyone interested in his own possible application. This statement is especially 
true for the photographic -detection mode. This note will attempt to elaborate on the 
method and its limitations. 
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The usefulness of wall -stabilized or cascade arcs as plasma light sources is well 
established. (For example, see ref. 1 (p. 205) and references therein.) Local thermo- 
dynamic equilibrium (LTE) and cylindrically symmetric geometry allow chordal observa- 
tions to be interpreted in terms of radial temperature profiles. Such an arc serves in 
this study as a convenient tool for studying the experimental aspects of Abel inversion. 
The working gas was atmospheric pressure argon and consequently, the temperature pro- 
files reported simply corroborate an existing bank of published data. (See refs. 2 to 5.) 

Several features distinguish this treatment from others. Instead of taking the Abel 
transform description as a starting point, as is usually done, it is developed here from 
basic considerations as the limit of a sum, and its solution is explicitly carried through 
to underscore boundary conditions entering into the problem. The two alternative forms 
to which numerical analysis is often applied are shown. 

The physical stops which limit spatial resolution are identified with due regard to 
their effects on depth of field and optical resolution. An approximate expression for the 
relative volume error incurred because of the use of converging optics is derived and is 
helpful in assessing alternative experimental schemes. 

Development of the Abel transform from a volume -weighted sum as is done here 
clarifies the optical thin assumption underlying the end result. The danger of relying too 
much on theoretical assurances of optical thinness is demonstrated by contrasting theo- 
retical and experimental estimates of the peak absorption for the 4300 Arl line. 


SYMBOLS 


A 


area, cm^; also mass, atomic mass units 


A u _^j transition probability, sec - * 


a 


line shape parameter, dimensionless, 


w col 

w Dop 


\J l°ge 2 


Bgtr brightness or steradiancy, used interchangeably with specific intensity. It 

is defined as the power per steradian per projected area of source in the 
direction of the detector, W(cm2-sr)"^. When subscripted with X it is 
a per spectral interval quantity having units Wfcm^-sr-A) - *. 


B functional form of steradiancy in terms of /3 having units of Bgj. r 

C empirical correction factor for intensity, dimensionless 
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optical density defined as minus the common logarithm of the transmission of 
a photographic negative or aperture diameter, mm 

Stark shift, angstroms 

energy level, cm"! 

focal length, cm 

oscillator strength, dimensionless 
dummy function 

apparatus function, angstrom" 1, or statistical weight, dimensionless 
slit height (that is, dimension perpendicular to direction of scan), mm 
sagitta, cm; also Planck's constant, 6.623 x 10"34 joule-sec 
specific intensity, W(cm2-sr)~* 

proportionality constant (deflection of densitometer pen per unit exposure 
on film), mm-erg"! (10^ ergs = 1 joule) 

linear absorption coefficient, cm"! 

line shape, angstrom"! 

lateral magnification, dimensionless 

number density, cm"3 

number of zones, dimensionless 

dummy variable, cm2 

radius at which arc ceases to emit detectable radiation, mm 


reflectivity, dimensionless 



r polar radial coordinate, cm 

r e classical radius of electron, 2.8179 X 10 - 13 cm 

S lens to arc distance, cm; also integrated absorption, dimensionless for 

wavelength-related case 

Sgig signal deflection, mm 

s geometric optics distance, cm 

T temperature, K 

t time, sec 

u,v dummy variables for integration by parts 

V volume, cm"3 

W slit width (that is, dimension along scan direction), mm 

w line width at half peak, angstroms 

tv dimensionless slit width 

Z partition function, dimensionless 

a scaling factor for polar radial coordinate, cm2 

0 scaling factor for Cartesian coordinate, cm2 

/3jj ijth element of Abel inversion matrix, dimensionless 

A depth of field, cm, or zone width for inversion, cm; also used as a small 

difference 

Ax,Ay,Az small increments of Cartesian coordinates x, y, and z, cm 
6 depth of focus, cm 
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6(x) Dirac delta function 

e volume emission coefficient, W(cm3-sr)"^; also minimum resolvable 

dimension on film, mm, or emissivity, dimensionless 

e functional form of the volume emission coefficient when written in terms 

of a, W(cm3-sr)~* 

v frequency, sec“l 

0 solid angle, sr 

r transmission, dimensionless 

9 angle defined in figure 11, radians 

£ Guant factor, dimensionless 

Subscripts: 

B brightness 

b back, that is, side furthest from detector 

c collecting 

col collision 

d densitometer 

Dop Doppler 

e electron 

eff effective 

f front, that is, side nearest detector 

1 lower level 
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o 


line center 


max maximum 

pyr pyrometer 

s slit 

t true 

u upper level 

X wavelength, a spectral quantity 

□ rectangular 


Superscripts: 

(arc) 

arc 

cont 

continuum 

L 

line 

meas 

measured 

nd 

neutral density 

0 

black body 

obs 

observed 

opt 

optics 

Pi 

plasma 

thin 

optically thin 

(W) 

tungsten 
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image space quantity 


+ ion 

DESCRIPTION OF THE ARC LIGHT SOURCE AND ITS OPERATION 

The light source is patterned after the design of Shumaker (ref. 6). The arc is 
maintained between two cylindrical tungsten electrodes which serve as cathode and anode 
and is constrained by a stack of annular, water-cooled segments, each having a central 
bore of 6 mm. The segments are separated from each other by the ring windows used 
for viewing and from the arc housing by dielectric spacers. The interelectrode spacing 
is fixed at 7.6 cm. 

The arc unit is pictured in figure 1(a) and a schematic labeling of the main compo- 
nents is shown in figure 1(b). The control center for the arc is shown in figure 2(a) and 
the flexibility of the gas porting system is suggested in figure 2(b). The option of 
exhausting into a manifold which is vacuum pumped allows the system to operate at total 
pressures down to 6000 newtons per meter^ (0.06 atmosphere). 

Electrical power is supplied by three arc welder units suitably filtered. For the 
source to be adequately described by a constant -temperature profile, good filtering is 
essential since the current density which controls it varies linearly with current in a 
wall -stabilized arc. Figure 3 shows the electrical schematic. Bus-bar arrangements 
allow machines to be connected with or without a water-cooled load resistor to yield a 
current range of 30 to 150 amperes in argon. The operating limit is set by total voltage 
available from the supplies. Considerably higher currents could be achieved by short- 
ening the arc column. Figure 4 depicts the arc characteristic curve. Total arc voltage 
was determined by measuring the field gradient between consecutive disks and multiplying 
by the electrode separation to guard against changes in the anode and cathode fall regions 
which can occur when increased current heats the electrode. (See ref. 7, p. 29.) 

To strike the arc, an end window is removed and voltage is applied to the elec- 
trodes. The operator deftly inserts and withdraws a carbon rod, touching the rear elec- 
trode, and subsequently draws the arc which is formed through to the near electrode. 

The arc usually spins around the cylindrical electrodes for several minutes and finally 
attaches to one spot and remains quiescent for hours thereafter. 

The arc runs best in pure argon as has been noted by previous investigators. The 
temperature dependence of its electrical and thermal conductivity is apparently well 
suited to the dimensions and current ranges of this experiment. In contrast, when helium 
or nitrogen is used in the arc, it becomes quite unstable and eventually attaches to one of 
the segments, and the run must be terminated to prevent cooling failure. Others have 
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(a) Side view of arc housing. 
Figure 1.- Stabilized arc. 







L- 71-5862 

(a) Current and gas flow controls. 

Figure 2.- Control center for the arc. 


10 


II III! 


in in 














so r 


Experimental point with estimated error. 



Arc current, amperes 

Figure 4.- Characteristic curve of the 7*6- cm arc run in atmospheric argon. 



experienced similar problems with hydrogen (ref. 8). To use this arc in such gases, it 
is usually advisable to isolate the test gas in the center section and use the argon to 
blanket the electrodes as suggested by Shumaker (ref. 6). The porting arrangement was 
designed with this procedure in mind. Because of the low flow rates used in this work, 
the arc is not well suited for introducing nongaseous contaminants. Attempts to introduce 
copper by placing a copper rod on one of the electrodes was unsuccessful. The blue- 
green radiation characteristic of copper could be seen next to the electrode, but evidently 
the copper plated out long before entering the test region since no trace of copper was 
detected spectroscopically. 

Since this configuration is a wall -stabilized configuration, it is natural that the arc 
will tend to bulge out into the viewing ports. The windows are 5 mm wide and the bulging 
is very noticeable, especially at high currents. The bulging often does not occur sym- 
metrically but will develop peculiar kinks because of some unknown combination of flow 
rates and field gradients. Any deviation from cylindrical symmetry negates the Abel 
inversion concept. Asymmetry can usually be detected a posteriori by folding the densi- 
tometered data about its center line to look for features on one side that do not appear on 
the other. To illustrate the effect, figure 5 was prepared from an isodensity trace of a 
plate taken in the slitless mode of a 3.4-meter Ebert spectrograph. The instrument 
aperture ratio is f/33; the instrument has a plate factor of 5 A/mm and 0.04 A resolution 
in the first order by use of a 590-groove/mm grating. (10 A = 1 nanometer.) The arc 
was viewed from the side so that the viewing port itself acted as the slit. The figure is 
interpreted in the following manner: The clear area represents background density on 
the plate. When the density increases by a predetermined AD, the plotting code changes 
in the order blank, dot, dash, blank, etc.; the order reverses if the increment is negative. 
Figure 5(a) shows a case which was judged visually to be symmetric and figure 5(b), one 
which was obviously distorted. 

Narrower windows, shorter arc column, and smaller bore should improve the per- 
formance of the facility with the harder -to-run gases and lessen the bulging problem in 
all cases. 


MEASUREMENT OF RADIAL EMISSION COEFFICIENT 

In an arc where the electric field defines a direction in space and many processes 
are governed by diffusion, there is a natural tendency for axial symmetry. This natural 
tendency is enhanced by the boundary conditions in the wall -stabilized or cascade arc. 

In this situation, it is assumed that the physical parameters characterizing the phenome- 
non (for example, emission coefficients, particle densities, temperatures, etc.) are func- 
tions of the radius r only. This dependence can be extracted from observations along 
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BORE DIAMETER -*) 




WINDOW WIDTH 



( LINE PLUS CONTINUUM ) 


(a) Symmetric conditions. 



(b) Asymmetric conditions. 

Figure 5-- Isodensity traces of arc column in vicinity of 4300 Arl line 

radiation on underlying continuum tends to define extent of window which 
is serving as a slit for this data. 
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chord length by mathematical techniques often referred to as inversions since they either 
involve finding inverse integral transforms or the inverse matrix of a linear 
transformation. 

In this section the Abel integral equation will be derived from basic energy and 
geometry considerations. The extent to which the requisite geometry can be approxi- 
mated in practice because of optical limitations is then discussed. 

Derivation of the Abel Equation 

When the collecting optics are highly collimated, the inversion problem can be 
stated and solved as an Abel integral transform. Consider the energy collected by a 
rectangular aperture viewing along a chord length as in figure 6. 


X 



(a) Cylindrical lamina showing second rectangular swath observed 
under ideal Abel conditions. 


x 
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(b) Axial view of cylindrical emitter. 
Figure 6.- Source geometry. 



Assume, subject to later analysis, that no convergence is introduced by the optics 
and, furthermore, that all volume elements are equally effective in contributing to the 
signal (optically thin implication). The aperture cannot detect the depth of the source 
and the intensity may be interpreted as emanating from a strip having a certain bright- 
ness or steradiancy. This approach is reasonable since absolute intensity is determined 
by comparison with the steradiancy of a tungsten filament. Let Bg tr x be the specific 
intensity or spectral steradiancy of the source. If e^(r) is the spectral emission coef- 
ficient, then for a slab of height Az as shown in figure 6 


B Str,xfcO Ax Az = ^ e x (r) Ax Az Ay 


( 1 ) 


-y f 


Since Ax Az is uniform and divides evenly from both sides, and reflection symmetry is 
used after proceeding to the limit of the sum, 


B Str,X« = 2 Iy = „° e X (r) ^ 


.1/2 


But y = (r^ - x^) 1 and, for fixed x, dy = r(r2 - x^) dr. Consequently, 


B Str,X 



e x (r) r dr 

(r 2 . x 2)V2 


( 2 ) 


Since the left-hand side is measured, one has an integral equation for the unknown quan- 
tity. This is a particular example of a class of integral equations known as Abel equa- 
tions. The kernel is singular because of the lower limit. The solution, although well 
known, will be outlined. 

Let a = r2 and jS=x2. If r is written in terms of a and substituted in e (r), 
a new functional form results which is designated by an overbar; hence, 


6 Or) -€(*) = 6- (r*) 


(3) 


and similarly, 


B Str,X^ “ B Str,X^ 


(4) 
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The integral equation is thus transformed to 


B Str,A 



e(a) da 
(a - /3)*/ 2 


Multiply both sides of this equation by d/3(/3 - p)”*/ 2 and integrate from /3= p to 
j8 = R 2 . Then 

P/3=R 2 Bgtr.xO) = d/3 f Q!=r2 6(a) da 

^/3=p (j8 - p) 1 / 2 ^/3=p (/3 - p)l/2 (a - 0) 1 / 2 



e (a) da 

(0 - p) 1/2 (a - /3)V2 


The integrand is some function of a and /3, for example, f(a,/3), and the order of inte- 
gration suggested is that /3 is held fixed and a varied from the line a = / 3 to 
a = R 2 as shown in figure 7(a). The same "area" is covered by adding strips at con- 
stant a as in figure 7 (b) with /3 going from p to a . Thus, 



(a) Summing vertical strips. 



Figure 7.- Strips used to define "area" of integration in the Abel derivation. 
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j^* 2 dg . f= R2 da C^ a ma) 

*^/3=p (/3 - p) 1 ' 2 ^a=p *^/3=p (j3 - p)V 2 ((j - jS) 1 / 2 


-r R2 m*»r Jfi — 

J C8=P J /3=p (j8 - p) 1/2 (0! - /3) 1 / 2 


The last integral is tabulated (ref. 9, no. 195.01) as 


2 tan- 1, /EL 

V« - / 3 


/3=oe 

i*=P 


= 2 tan"l(+°°) - 2 tan" 1(0) = 77 


By use of this result, the equation reads 

r /3=R 2 Bg t x (/3) pQ!=R 2 

- Str ^ 1/9 ^=77 e(a)da 

J /£=P (jS - p) 1 / 2 


'Q!=p 


Differentiate both sides with respect to p to obtain 
d C ^ r2 B Str,X^ Hp - f v 

TT \ 7 \i J 0 d/3 = "776 (p) 

d P ^ /3=P (0 - p) l/ 2 



(5) 


Let p = r 2 and recover the initial variables to obtain 


e(r) = 


1 d 
27 tt dr 


r- 


•x=R 


L 


x=r 


B Str.X (x) 
(x 2 - r 2 ) 1/2 


x dx 


(5a) 


This form is used by Barr (ref. 10) as a starting point for his numerical evaluation since 
smoothing can be done prior to numerical differentiation. The more common form is 
obtained by integrating equation (5) by parts before differentiating. Let 


u=B Str,X 

dv = (0 - p)" 1 / 2 d/3 
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Then 


dU = %^ 


v=2(0-p) 1/2 


and 


e -<p)=-±f 

Ti dp 


1 _d_ 
* dp 


2B Str,A^> ~ P) 


1/2 


/3*R^ 

;3=p 


-if'** o- P )V2^A d(? 

J^=p dp 


dB c 


2 B str,x< R2 > < r2 ' e > 1/2 - 2 I^ p R ((3 - »> 1/2 d(3 


dB c 


But B™. , (R 2 ) = Be- . (R) = 0 since R is defined as the point at which the signal is 

OCI*jA btl* j A 

zero. Then 




1 r^ R 1 ^Str,A 

17 ^P=p (/3 - p) 1 / 2 

Again letting p = r 2 and using original variables yields 


dB c 


e(r 2 ) s e(r) = 


j pX=R 

77 *^x=r 


dB 


(x 2 - r 2 ) 1/2 


Str^dx 


dx 


( 6 ) 


Although many papers have appeared on the numerical evaluation of this integral, 
only two will be cited. Bockasten (ref. 11) evaluates equation (6) directly in a rigorous 
fashion. Since experimental data are numerically differentiated in this method, fine 
zoning and low-noise data are required. Barr’s method (ref. 10) of working with equa- 
tion (5a) and smoothing the data prior to differentiation is used in this work. Barr com- 
pares the two approaches in his article. 
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Optical Considerations for Abel Inversion 

To utilize these analytical results, the optics are arranged to project the small 
rectangular field stop of the system into the emitting source so that energy will be 
accepted from small swaths (see fig. 6(a)) along successive chord lengths of the source. 
Apertures in the system should be small enough to provide as many nonoverlapping 
swaths as possible since one is essentially representing a continuous function by a series 
of discrete points. The stops and optics quality must insure sharp focus for all points 
along the diameter of the source; otherwise, the volume elements would be unequally 
weighted contrary to the underlying assumption. In practice, stopping the system down 
to achieve the requisite depth of field and data point density is impeded by the onset of 
diffraction and/or poor signal-noise ratio. 

Most published information relating to radial profile measurements deals with 
photoelectric output and the scanning slit aperture is located at the stigmatic plane of the 
spectrograph (entrance slit). When the photographic mode is used, the densitometer scan- 
ning spot as imaged into the source ultimately becomes the field stop. Consideration of 
figure 8 should clarify this point. The densitometer scanning beam with dimensions 
Wjj by intercepts an area of the processed plate. The transmission of the processed 
plate is uniquely related, through sensitometry, to the illuminance of the plate during its 
exposure. By neglecting the spread function of the emulsion in comparison with scanning 
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slit dimensions, the latter can be thought of as imaged by the spectrograph and inversion 
optics into the volume emitter where it has a cross section by H^j. 

To assess the stigmatic quality of the spectrograph, a reticle that is graduated 
finer than the desired zone spacing is placed at the entrance slit. Its image is viewed 
in the film plane and examined for distortion. This examination is made either in white 
light by using zeroth order or in monochromatic light (for example, Hg green line) by 
using the dispersion mode. 

The optical layout is sketched in figure 9(a) and pictured in figure 9(b). The 
"T" arrangement allows for absolute intensity calibration with good repeatability and 
minimal effort. The arc is mounted on a rotatable base so that it can be viewed axially 
or side-on. The optical train consists of a single lens L and aperture of diameter D. 
This aperture controls the energy throughput (aperture stop), degree of parallelism, and 
depth of field. The depth of field dependence on D will now be considered, the analysis 
of the convergence aspects being deferred to the next section. 

Let e be the minimum resolvable dimension on the film or perhaps, to be safe, 
some multiple of it, and S' the lens-to-slit plane distance. Then with m, the linear 
magnification of the lens system, the following relationship for the depth of field A is 
derived in appendix A: 


A 


1 2S’e 
m 2 D 


In practice, the source is often very small and no significant demagnification can be tol- 
erated. The quantity S' is reasonably determined by spatial requirements of the labo- 
ratory and e does not vary significantly for available spectroscopic plates. By elimi- 
nation then, the depth of field is mainly regulated by the lens of diameter D. If 2R is 
the diameter of the arc column, then 


1 2S’e 
m 2 D 


s 2R 


or 


m 2 R 

The diameter D must not be reduced to the point where diffraction deteriorates the 
image of the arc column on the slit. For arcs, it is feasible to look for diffraction effects 
by placing a white card at the entrance plane. For nonvisible light or weak sources, one 
might insure that AS' < (ref. 12). 
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3.4 - METER 
EBERT 

SPECTROGRAPH 


(a) Optical layout. 



L-71-5864 

(b) Partial view of the experimental layout with the arc in the end- on position. 


Figure 9*- Experimental setup. 
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Estimate of Nonparallelism 


In figure 10, a typical volume of intersection is sketched with convergence present. 
It should be contrasted with figure 6(a). An expression for the effective slit width W e ff 
will now be derived. The effective width on the side of the source nearest the spectro- 
graph is called the front effective width and would be subscripted with an f to contrast 


ARC 

AXIS 

Z 



Figure 10.- Sketch of prismoid of intersection for actual case having finite f numbers. 

it with the back effective width W eff b and usually the subscript eff will be suppressed 
to condense notation. The width W b is always larger and will determine the distance 
between observation points and hence the number of zones possible with a given geometry. 
The result will be useful in estimating the relative volume error made when the Abel 
requirements are compromised. 

In figure 11(a) the extreme rays are sketched for one dimension of the scanning slit 
when it is centrally positioned. The relative dimensions are distorted to emphasize con- 
vergence. Clearly, a conservative estimate for W b is the chord whose half length 
is ab. In figure 11(b) construction lines and symbols used to relate ab to known quan- 
tities are provided. 
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FIELD STOP 


SOURCE 


¥ 

APERTURED LENS 

(a) Extreme rays for swath 1. 



(b) Lines used in derivation of W e ff t 
Figure 11.- Derivation of 


As a result, 


W e ff,b = W' + 2y 


Considering triangles 123 and 345 (fig. 11(b)), 


D . W 



or 


y = 


D W1 
2 2 
S’ 


x 


By noting that 


x = R cos 9 = R(1 - sin2$) 


1/2 


9 


w b = W' + P + W * R(1 - 

s 

( w \ ^ 

— Jh = — and 
2R/ 4n 2 

Even for a relatively coarse five zone case, the radical term is approximately 1 and need 
not be carried 


W v 


w , R(D + W’) 
S’ 


The quantity Wf is calculated in a similar fashion and introduces a minus sign in the 
parentheses and S' must become infinite to recover the Abel geometry. In practice 
S' is made large and D small, subject to the constraints listed in the preceding section. 
It is desirable to have an indication as to the effects that the finite optics will have on 
interpretations of emission coefficients. By using W b as zone width and interpreting 
the data with Abel equations, the source volume collected is overestimated. Consequently, 
the energy per unit volume which is calculated will be lower than the true value. Further- 
more, the error should be largest in the central hot zones where the volume is smallest. 

A precise calculation of the volume omissions involve triple integrals with limits set by 
the intersections of prismoids and cylinders as suggested in figure 10. 


Since writing down the expression for the volume integrals is not particularly 
instructive, the following more tractable estimate is suggested. A volume element 
bounded from below by dA and from above by f(x,y) is f(x,y) dA, and for finite volume 


V = 


If 


dA. 


If f is a slowly varying function of x and y as it would be on one 


of the faces of the prismoids, then f could be taken from the integral as approximately 
constant. In this spirit, consider figure 12(a) where a swath is depicted in projection and 
an enlarged view of the zone 1 in figure 12(b). The term W' is assumed to be centered 
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(a) Figure used to estimate relative error in volume measurement 

due to finite optics. 



Figure 12.- Estimation of relative error. 


in swath 1 for simplicity. It is apparent that the area omitted will be less than or equal 
to the hatched area Aj + A 2 . Then, reasoning as above, the relative volume error is 
taken as 

AV _ 2AA < 2 ( A 1 + A 2) 
v 7rW b 2 " 7rW b 2 
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From plane geometry (ref. 13, p. 12), recall that the area of a segment bounded by 
a chord a distance R - h from the center of a circle having radius R is 

r 2 cos -l^_li^ _ (r _ h)(2Rh - h 2 ) 1 / 2 . Since hj = (W b + W')/2, then 


A 1 = W b‘ 



+ ^3 2 + V 

2 1/4 b 2 


Similarly, h 2 = (W^ - W’)/ 2 giving 

0 /W K + W’\ W, +W' f 

cos ‘\-lwg-) - -V-] V; 


W 


,2 


3 „ T 2 W b W ' W 2 
^ W b 2 4“ 


W T 

Introduce these expressions in the relation AV/V and define w s — - which is the 

w b 

reciprocal of the effective slit dimension in units of its geometric image as determined 
from optics and obtain 


1 + 2 
V 7T 


cos ‘HH^) * cos 'K^) 
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(1 - zo)\J (3 - u/)(l + to) - (1 + ^)/(3 + zu){l - to) 


By using the approximate expression for one obtains 


w 


.(i.S-it'f 

V w s’ s’/ 


The dimensionless ratios R/S' and D/W* quite clearly control the relative error. 

Note that if S' becomes infinite, w becomes 1 and the relative volume error is zero 
as it should be for the Abel case. If a rectangular scanning aperture is used, the expres- 
sion should be tested for both slit dimensions. To guard against convergence problems, 
R/S' and D/W should be chosen to make AV/V less than 5 or 10 percent. Fig- 
ure 13 shows the dependence of AV/V on both D/W' and R/S'. For the data reported 
in this note D/W' was 1.6 mm/1.6 mm or 1, and R/S' was 3 mm/1000 mm or 0.003. 

It should be emphasized that this expression is approximate and that the curves should 
not be used to correct measured values. 
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OPTICAL DEPTH CONSIDERATIONS 


The solution of the radiative transfer equation for slab geometry in a collision- 
dominated gas having a Planckian source function (based on ref. 7, eq. (69.9)) is given by 

- exp(-k x x)J 

where Bg tr X (T) is the specific intensity of a black body at temperature T, x the 
geometric depth into the slab, and k x is the spectral linear absorption coefficient. It 
is tacitly assumed when inverting the data that each photon emitted in the source escapes 
from it without reabsorption. Mathematically, this relationship holds if the optical 
depth k x x is small compared with 1 so that two terms in the Maclaurin expansion of the 
exponential suffice to represent it, namely, 

^ in “B^ tr>x (T)[l-(l-k x x)] 

for then the signal increases in direct proportion to the depth of source viewed. 

Theory 

The line is fully characterized by its absorption coefficient, which it is convenient 
to write as 



k x = Z(\ - A 0 )S 

Thus, spectral detail as carried by the normalized line shape X, that is, 
f°° X(X - A 0 ) d(A - Ao) = 1 

_ oo 

is explicitly separated from the strength of absorption determined by the integrated 
absorption S. An explicit expression for S when written in terms of wavelength is 

S = A 0 2 7rr e N 1 f 1 _ u 

The specific intensity of each line "grows” with optical depth in its own peculiar 
way which must be determined empirically and results in its curve of growth (ref. 14, 
p. 64). The region of this curve which is linear in x identifies the optically thin domain. 
In practice, it is often sufficient to calculate k x by using an approximate model to 
insure that the truncated series expansion is valid. A popular expression is the Voigt 


30 



I 


form (ref. 15, eq. (3-34)) which is the convolution of dispersion and Gaussian shapes 
parameterized by 

a = ^““( Iog e 2 ) 1/2 
Dop 

where w cQ ^ and w j) Q p are, respectively, the full widths at half maximum for the col- 
lision and Doppler broadened components. The Gaussian profile is recovered from the 
Voigt formula when a is zero, and the dispersion profile when a is infinite. Although 
the Voigt formula incorporates the overall general aspects of line broadening, its approx- 
imate nature should be stressed. Its symmetric form suggests that only those collisions 
which terminate the emitting wave train are included in the collision broadening, whereas 
observed Stark profiles display a strong asymmetry and shift attesting to the vast number 
of the collisions which only alter the phase of the wave train. To show its usefulness, a 
rough estimate of the peak absorption expected for a line in the plasma will be given. 

To be specific, consider the 4300 Arl line as radiated from a 10 000 K argon plasma 
at a total pressure of 1 atmosphere. Penner (based on ref. 15, eq. (4-37)) shows that 

k X Q = exp(a 2 )[l - erf (a)] 

The composition data of Popenoe and Shumaker (ref. 16) were used to obtain table I. If X Q 
is the line center in angstroms, T, the absolute temperature, and A, the mass in atomic 
mass units, then the Doppler half intensity width is given as (ref. 17, p. 67) 


w Dop = 7 - 11 x 10_7 /i X o’ angstroms 

= 7.11 x lO-y^OO 43 qq _ o.048 angstrom 

The collision broadening is mainly electronic Stark broadening which is obtained from 
Griem (ref. 1, p. 492), and by referring to table I for the electron density 

2w(10 000 K)N 

w col - t - 10 16 ' angStr0mS 

= 2(0.106)(l.511xlQl6) _ ^ ^ 

1 X 10 16 
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The a parameter is now easily evaluated as 5.6. Since the absorption oscillator 
strength is related to the transition probability through (ref. 18) 

*l-a= 

with X in angstroms, then using 0.411 x 10® ± 11 percent sec - l for A u _^ (ref. 11) 
f£00 Arl = ! 4992 x 10“ 16 (4300) 2 1(0.411 x 10®) = 0.0019 


The integrated absorption then becomes 

S = (4.3 x 10"5) 2 (3.14159)(2.8179 x 10-13)(2.930 x 10l2)(0.0019) 
= 9.1 x 10" 12 (dimensionless) 


erf (5. 6 

= 0.0018 cm" 1 


giving the peak linear absorption coefficient estimate of 


k 4300 Arl . 2(9.1 X 10 ^, 1 ^ CgSl r . 

X ° 0.048 x 10-8 V * L J L 


Since x max is on the order of the diameter of the column for the inversion work, 

2Rk^ o « 0.001 « 1 and the optical thin assumption seems to be well justified. 

Experiment 

The arc configuration allows for empirical verification of the optical depth estimate 
and, more significantly, for a direct evaluation when line-shape information is unavailable. 
The arc chamber is rotated so that it can be viewed axially. By looking at a small cross 
section, isothermal zones of plasma are viewed. A concave mirror is placed behind the 
column at a distance 2F from the arc midplane where F is the focal length of the 
mirror and in such a way that, when uncovered, it directs the arc radiation incident upon 
it back through the system and onto the entrance slit. The mirror is alined by directing 
a He-Ne laser beam through the film plane of the spectrograph. It is reflected from the 
internal optics and out the entrance slit, down the optical axis of the system, and through 
the arc chamber to the concave mirror. The mirror is adjusted to cast the beam back 
through the entrance slit. Since the optical lever arms are on the order of several 
meters, it is felt that this arrangement insures very good alinement. From consideration 
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of the sketch of figure 14 with the end windows having transmission r, the mirror having 
reflectivity It, and tP* = e - ^ representing the plasma transmission, it is evident that 
the ratio of double-path to single-path signals should be 1 + t^tP*, all quantities being 
spectral. Figure 14 reproduces an actual data trace. At selected points, the single- 
path signal is multiplied by 1 + Xt% and is plotted with error bars estimated from 
repeatability tests. Since the reflectivity and transmission values are known, the ratios 
of double -path to single -path signals from figure 14 allows the determination of rP* 
which can be solved for k^. Spectral resolution was adequate enough to allow measure- 
ment of peak line values, and k^ Q for 4300 Arl was measured as 0.14 per cm and con- 
firms the optical thin analysis. The discrepancies in k^ Q however, underscore the 
caution of relying over much on estimates using approximate line shapes. 

The usefulness of the experimental arrangement will be most evident when constit- 
uents are tested which have a large spread in lower energy levels and whose spectrum 
should therefore display a wider variation of spectral absorptivity than is afforded by 
argon. 


TEMPERATURE MEASUREMENTS 


Absolute Line Intensity Method 

The temperature measurement rests ultimately on a comparison of the brightness 
or steradiancy of an atomic line with that of a tungsten filament, the latter calculable 
from observing its brightness temperature. Whereas the total power output of the 
line 1^ is understandable from energy principles alone, a calculation of its spectral 
steradiancy requires a knowledge of its line shape. It is not immediately obvious how 
this strongly wavelength-dependent term is to be measured. 

Consider the problem of measuring the spectral specific intensity 1^. It is useful 
to distinguish further the true quantity as would be detected with infinite resolution from 
the measured quantity under finite resolution. This distinction is made with appended 
superscripts. A spectrograph distributes a monochromatic beam of light of wave- 
length A 0 over a range of wavelengths whose functional dependence is designated by 
g(A - A 0 ) and is called the apparatus function. The true and measured distributions are 
related through a convolution; namely, 


T L,meas 

L X 


= f I^,’ trUe g(A - A') dA’ 

_ OO A 


The apparatus function is determined by scanning a very narrow atomic line as would be 
emitted, for example, from a low-pressure microwave discharge. Then 
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jL,true _ (constant) |s(X’ - X 0 jj 


and 

jL,meas _ (c 0 nstant)|g(X - X Q )J 

The function g is usually normalized so that its effect will be to redistribute 
while conserving energy. If the entrance slit of a spectrograph is very wide, a rectangu- 
lar apparatus function g Q obtains 

S D ( x - x o) = ^ (| x - x o| s t) 

g D (x-x 0 ) = ° (l x - x °|>|) 

where W is the spectral slit width obtained by multiplying the slit dimension by the 
reciprocal linear dispersion of the instrument. To this approximation 

I L,meas = C°° (x _ v)I L,true dx , 

A J-OO D A 

= ±. r w/2 i^ true dx- 
w J-w/2 X 

and if the extent of the line is much less than W, little accuracy is lost by extending the 
limits to infinity to yield 

T L,meas _ 1 f T L,true _ I ’ 

: x ~wj. co y W 

To summarize, whenever the spectral slit width can be made much wider than the 
true line width to justify the redefined limits, the total line output follows by multiplying 
the spectral quantity at any point within the signal by W. Unfortunately, for the strongly 
broadened lines in plasmas, this must be done with some reservation. In figure 15 a 
trace of a mercury line suitable for apparatus function determination is shown along with 
an argon line as it appears in an arc plasma. Both traces are taken with the same spec- 
trograph settings and differ only in that the arc line is not narrow with respect to the 
spectral slit width - it is, in fact, comparable with it. 

To solve this problem, a rectangular apparatus function is generated by using wide 
slits, and the total intensity derived on this basis is empirically corrected. This 
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procedure will now be described. Separating the defining integral for total intensity I L 
into parts gives 

jL _ C°° jL.meas^ _ P'W/^L.meas^ P w /2 ji^meas 
J»oo * J_oo * J-w/2 * 

The first and third integrals on the right-hand side are designated by AA. It is conven- 
ient to center the integrals at A 0 . In figure 16 the different contributions are identified 
by hatching. Thus, 

I L = f W/2 lJ J » meas dA + AA 
J -W/2 A 

and by making the rectangular apparatus function assumption, 

I L = W . j^ J > meas _ A A n + AA 
A o u 

= W • I^meas + M = ^ + a CJ L 

where W • I^ meas has been abbreviated to Ip, and the last equality defines the correc- 
tion factor C. Figure 16 also identifies AI in terms of areas under observed traces. 

To reiterate, W, the spectral slit width, is obtained by tracing a narrow line from a low- 
pressure discharge. The arc line is read up by using a rectangular apparatus function 
and corrected by using the empirical C determined by planimetry. The empirical C 
must be redetermined for each line and whenever conditions are expected to change the 
line shape. At high currents, the correction can amount to several hundred degrees. 

The procedural aspects of "inverting” photographically have been commented upon 
with regard to figure 8. It is now appropriate to enlarge that discussion to indicate how 
the absolute intensity information is obtained. The cylindrical arc column is imaged on 
and perpendicular to the entrance slit having dimensions W s by H g . The instrument 
introduces no demagnification and is sufficiently stigmatic; therefore, the image on the 
plate in the vicinity of a spectral line will roughly have dimensions W s by H s . The edge 
will be diffused along the dispersion dimension because of the nonmonochromatic nature 
of the line which was corrected for above and along the height dimension because the 
radiance of the arc decreases as chords closer to the outer edge are viewed. To simplify 
the analysis, figure 17 portrays the density of developed grains by a contour map. Also 
shown are the dimensions and scan direction of the densitometer slit W d by H d . The 
dimension H d is taken small compared with W s so that the trace will be related to 


dA + r iL,mea Sdx 
J W /2 A 
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Figure 17.- Contours of developed line showing photographic inversion technique. 


the spectral intensity at A 0 . The dimension may be adjusted within the limits set 

by the desired inversion resolution and adequate signal-noise ratio. A spectral sensi- 
tometry strip is recorded on each plate to relate the densitometer signals to relative 
intensity. Care is taken to maintain chemical purity and temperature control during 
processing and a nitrogen burst system is used to agitate the developer as a safeguard 
against the Eberhard effect. The arc is much more intense than the tungsten filament to 
which it is compared and its exposure time orders of magnitude shorter. To guard 
against reciprocity failure, a neutral density filter having transmission r nc * is intro- 
duced into the optic train for the arc exposure, its magnitude being chosen to make expo- 
sure times equal. 

Reconsider figure 9(a) to examine the irradiance of the film. The aperture or 
energy -limiting stop for the system is the circular mask on lens L. The numerical 
aperture of the spectrograph is small enough to permit all energy passing through the 
slit to be gathered by the collimating mirror. First, consider the signal obtained when 
the pivot mirror is rotated to accept energy from the tungsten lamp whose specific inten- 
sity is ijW (Tt) where T^ is the true temperature. The solid angle of acceptance as 
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determined by the aperture will be designated S2 C . The densitometer scanning aperture 
serves as the field stop and consequently determines how much of the source area is 
viewed. With lateral magnification m defined as before, the source area viewed is 
seen to be Let be the proportionality constant relating the exposure of 

the film to deflection of the recording pen. It will be a nonlinear function involving the 
film characteristic, electronics, and so forth, but is known from sensitometry. Then, the 

signal S® is 


s<w> = 

sig 

■*(■3?) 

opt 

T X 

W — 
s dx 

,(W) 

W d H d 

m^ 

«c 

lf>(T t ) 

Deflection 

Deflection 

per 

unit exposure 

Optics 

transmission 

Wavelength 

interval 

Exposure 

time 

Source 

area 

Solid 

angle 

Specific 

intensity 


It is possible to write r opt explicitly as the product of all lens transmissions and 
reflectivities of the system including grating efficiency. This is impractical and unnec- 
essary since, for a T system, the term cancels, as shown subsequently. With the mir- 
ror accepting arc radiation and with the addition of a neutral density filter 

4ig 0) * K x(4tg c) ) TOpt ^ dw B s t<arc) “4r «y| arc) c« 


As previously mentioned, 


T n d is adjusted to make t^^ = t( arc ^. Then 


T (arc) _ 
L X 


(arc) 
sig 



«x(»g?) igV> 

K x( s ^ c) ) 


The specific intensity of the filament is calculated from an observation of its 
brightness temperature by using a brightness pyrometer having effective wavelength 
Apyr = 0.635 micrometer. The pyrometer is placed for convenience as in figure 9. 
Then the brightness temperature observed depends on the reflectivity of the mirror % 
and Tg is obtained by iterating 

B Str,Apy r ( T B bS ) = ^A pyr B Str,A pyr ( T B) 
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The true temperature of the lamp is found by using emissivities from Vos (ref. 19) by 
iterating 


■p>° 

Str,Xpy r 


(T t ) = ejM 
A pyr 


T)p 

Str,Xpy r 


( t b) 


and, finally, the specific intensity at the wavelength of interest by 


(W) = e (W) R o 


Str,A 


(T t ) 


This reduction procedure must be repeated for each data point required. Further- 
more, in an arc plasma, the line rides on a continuum which can equal it in magnitude. 
Fortunately, the continuum is nearly constant and its component can be subtracted by 
using a trace adjacent to the line in question. Since the line intensity result depends on 
the differences of comparably large numbers and must (unlike the continuum) be corrected 
for line shape in addition, it loses some of its attractiveness as a temperature indicator. 
This problem is considered again in the next section. In figure 18 an actual block of raw 
data is shown with traces and zone boundaries labeled. Barr (ref. 10) has reduced the 
smoothing and Abel inversion mathematics to a matrix multiplication operation. The 
measured steradiancy Bg^ r ^(x) is approximated by a discrete data vector whose com- 
ponents are made up of the values taken by Bg^ ^ for equispaced intervals along x. 
Small segments of this data vector are then fitted with polynomials, substituted into equa- 
tion (5a) and subsequently integrated term by term. Let A be the interval spacing and 
Ij = Bg tr ^(iA); then Barr incorporates all of the curve fitting, smoothing, and integration 
into a matrix such that 

e i (r) = Vl 
j 

where e^r) = e (iA), that is, the emission coefficient for the ith radial zone. For a line 
having transition probability the power per unit volume per steradian is given by 


* - r, NuAu-l^lu 

where i/j is the frequency and h is Planck's constant. By using the composition data 
of table I and the transition probability of reference 16, the emission coefficient has been 
calculated and is plotted in figure 19. The radial temperature profiles derived by use of 
4300 Arl radiation are shown in figure 20. 
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Emission coefficient, |i watts/ cm' -sr 



Temperature, K 

Figure 19*- Variation of emission coefficient with temperature for 4300 Arl. 
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files for the 30- , 90 - , and 150- ampere runs obtained from 
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All the calculations are done with a machine program so that only certain constants 
related to the readup need to be entered along with direct deflections from the densitom- 
eter traces. In this context, only one comment will be made concerning the analytic form 
chosen for the film characteristic curve. The conventional plot of density against the 
logarithm of exposure results. in a more or less S-shaped plot referred to as the 
H and D curve after Hurter and Driffield. This form is an inconvenient one to represent 
analytically and a recent paper by de Vaucouleurs (ref. 20) redirected attention to a much 
earlier work by Baker (ref. 21) showing that the "toe" of the curve is considerably 
straightened by plotting log 10 (lO densit y - l) instead of density along the ordinate. For 
densities of 1 or greater, the curve practically reduces to the original H and D curve, but 
on the low end considerable linearization occurs. The sensitometry data represented in 
this way is linear least-squares fitted. The standard output lists the deviations from the 
fit and thus warns the user of anomalies in the data. 

Absolute Continuum Method 

Continuum radiation results from free-bound and free-free transitions in the arc 
plasma. The exact calculation of continuum radiation is only possible for hydrogenic 
systems. Chapter 5 of reference 1 provides a survey of the working rules for its calcu- 
lation. Simplified theory predicts the shape of the spectrum very well but leaves the 
magnitude uncertain by a proportionality factor involving radial overlap integrals called 
the Gaunt factor, |(X). This factor is often taken as a variable and adjusted empirically 
as is the case for argon. Corliss and Shumaker (ref. 22) provide the expression 

e cont . 9 ,8 x 10 ~ 21 *< X > V uwatts _ 

rp 1/2^2 g cm 3 -sr-A 

with X in angstroms, Z + the ion partition function, and £(4315 A) = 2.3. This expres- 
sion is used to calculate the last column in table I which is plotted in figure 21. The tem- 
perature profiles derived from continuum intensity measurements are shown in figure 22. 
The continuum measurements have the advantage over line measurements in that no line 
shape or background corrections need to be made. 

Electron Density Method 

The electron density N e is a unique function of temperature in an equilibrium 
calculation. This quantity may be used along with table I to ascertain the temperature. 
The shift and broadening of spectral lines are functions of the local electron density. In 
figure 23 is shown the argon spectrum at 30 and 150 amperes bracketed by an iron (Fe) 
reference. The shift and broadening of the high current spectrum is readily apparent as 
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Figure 22 .- Temper at ure profiles obtained from absolute continuum 
intensity at 4^15 angstroms. 
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L- 72-124 

Figure 23-- 30- and 150-ampere runs to illustrate enhanced continuum, line shift, 
and. line breadth at high electron density. 

p 

is the increased continuum contribution which goes as N e . The high-amperage case 
has been exposed slightly longer to enhance the asymmetry. Although some approximate 
calculations of shift d and full width at half peak w are tabulated (for example, ref. 1), 
again the safest procedure is to use empirical results obtained by performing experiments 
with traces of hydrogen added whose line shapes are very reliable and serve to determine 
the local electron density. Popenoe and Shumaker (ref. 16) have used this approach to 
determine constants pertinent to 4300 Arl and find 

w = 0.2 x 10" N e (A) 


and 


d = 0.103 x 10" 16 N e (A) 
with N e in cm"3. 

To avoid problems with unfolding line profiles, the arc is viewed axially and only 
the central zone is measured. The slits are set to diffraction fill the collimator. This 
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setting results in an essentially Gaussian apparatus function having a full width of 
0.016 angstrom. The narrowest full width at half peak measured was on the order of 
0.2 angstrom which is more than 4 times the spectral slit width of the spectrograph 
having a 10 -2 -mm (10-micron) entrance slit. Kostkowski and Bass (ref. 23, fig. 2) plot 
the relative width error encountered because of apparatus broadening as a function of the 
ratio of spectral width to true width. From this curve, the largest expected error encoun- 
tered in this experiment should be less than 2 percent and is neglected. 

Axial temperatures measured by use of Stark results are compared with those 
derived from other measurements in figure 24. The lines are very broad and the peaks 
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Figure 24,- Comparison of axial temperatures ascertained by different techniques. 


are difficult to determine; thus, the shifts are difficult to measure accurately. The width 
measurement is a little easier to do. Agreement seems adequate in view of the fact that 
formulas provided by Griem (ref. 1) for the same line differ from these values by 
20 or 30 percent. 
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CONCLUSIONS 


The Abel integral relationship for inverting "side on" spectroscopic data is seen to 
follow from assuming cylindrical symmetry, optical thinness, and nonconverging optics. 
Its restriction to sources which have a finite extent is seen to follow from the boundary 
conditions imposed on the explicit solution of the integral equation for the transform part- 
ner. From the analysis the following conclusions are made: 

1. The dimensions of the scanning aperture on the densitometer and the aperture 
stop on the lens which images the source onto the spectrograph slit combine to determine 
the spatial resolution obtainable with the Abel inversion. 

2. There will always be some underestimation of the emission coefficient derived 
from inverted data because of converging optics but this value can be kept small by 
appropriate choice of stops and optical component spacings. 

3. Approximate theoretical estimates of optical depth are risky and it is best to rely 
on empirical determinations if the source geometry allows for it. 

4. Absolute intensities for highly broadened lines require correction due to back- 
ground and line shape; thus, absolute continuum measurements are made more attractive 
for temperature determination when empirical Gaunt factors exist. 

5. Stark width and shift temperature measurements operate indirectly through the 
electron density, are difficult to acquire especially in the nonhomogeneous case, and 
except for hydrogenic systems are hard to interpret because of discrepancies in the 
literature. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 9, 1972. 
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APPENDIX A 

DEPTH OF FIELD FOR A SINGLE LENS 

Although most people have an intuitive feeling for depth of field, it is surprisingly 
difficult to come by its discussion in optics books and even more so to find a working 
formula for it. Consider figure 25 which depicts an ideal thin lens L having focal 


L,F 



Figure 25.- Single lens depth of field diagram. 


length F and aperture D. This lens will image an axial point at S into S’. The 
depth of focus concept is introduced as follows: If the plane of observation were moved 
by 6/2 toward or away from the lens, the spot of diameter e would still be an accept- 
able image of the point. As with any discussion of resolution, the criterion of accept- 
ability has an element of arbitrariness, but, in practice, it could be the minimum dimen- 
sion resolvable by a particular film as stated in the manufacturer's specifications. 

By using similar triangles 125 and 134 (fig. 25), 

D/2 _ e/2 
S' 6/2 

or 

6 = (Depth of focus) (Al) 

All axial points between Sj and Sg conjugate to s'j and Sg, respectively, will satisfy 
the criterion of acceptability and, in fact, s^ - Sg = A will be termed the depth of field. 
Since s'j = S' - & and Sg = S' + the Gaussian equations pertinent to the thin lens 
read 
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APPENDIX A - Continued 



Equating the left-hand sides of these equations yields 



and rearranging yields 


(A2a) 


(A2b) 


s 2 



1 



By using the least common denominators and the definition of depth of field A, 


A = 


6 


Ol2 6 2 

s "T 


S 1 S 2 


(A3a) 


A = 



s l s 2 


In practice, 


(e/D) 2 « 1 


and can be neglected. Then, 


A 


5 

S ’ 2 


s l s 2 


If Sj « s 2 ~ S, 


A 


6 2S'e 
m 2 m 2 D 


(A3b) 


(A3c) 


(A3d) 


gt 

where m is the linear magnification, - ^ , and its square is usually referred to as Ion- 

s 

gitudinal magnification. 
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APPENDIX A - Concluded 


Although this final form is concise and useful, it should be remembered that it is 
essentially a differential approximation and if A becomes, for example, 10 percent 
of S, it is advisable to use the relations (A3c) and (A3b) with and s 2 substituted 
from equation (A2). 
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TABLE I.- TEMPERATURE DEPENDENCE OF QUANTITIES ENTERING 
INTO THE SPECTROSCOPIC MEASUREMENT 

Composition data and spectroscopic constants from reference 16. 

4300 Arl: A = 4.11 x 10 5 sec-1 (±11%) 
g u = 5; E u = 116999.4 cm~l; 
gj = 3; Ej *= 93750.6 cm“l. 


Temperature, 

N e , 

N P 

N u , 

o 

e L 

c i 

H watts 

e cont 
p watts 

K 

cm"3 

cm -3 

cm"3 


cm^-sr 

(a) 

cm3-sr-A 

(b) 

7 500 

6.684 x 10 14 

4.540 x 10 10 

8.755 x 10 8 


1.324 x 10 1 

1.100 x 10° 

8 000 

1.453 x 10 15 

1.307 x 10 11 

3.331 x 10 9 

0.5 

5.037 x 10 1 

5.100 x 10° 

8 500 

2.898 x 10 15 

3.305 x 10 11 

1.077 x 10 10 


1.628 x 10 2 

1.980 x 10 1 

9 000 

5.352 x 10 15 

7.490 x 10 11 

3.037 x 10 10 


4.592 x 10 2 

6.520 x 10 1 

9 500 

9.254 x 10 15 

1.544 x 10 12 

7.614 x 10 10 


1.151 x 10 3 

1.892 x 10 2 

10 000 

1.511 x 10 16 

2.930 x 10 12 

1.723 x 10 11 

5.6 

2.605 x 10 3 

4.901 x 10 2 

10 500 

2.342 x 10 16 

5.163 x 10 12 

3.560 x 10 11 


5.383 x 10 3 

1.146 x 10 3 

11 000 

3.466 x 10 16 

8.488 x 10 12 

6.764 X 10 11 


1.023 x 10 4 

2.446 x 10 3 

11 500 

4.910 x 10 16 

1.307 x 10 13 

1.189 x 10 12 


1.798 x 10 4 

4.789 x 10 3 

12 000 

6.672 x 10 16 

1.887 x 10 13 

1.937 x 10 12 

23.7 

2.929 x 10 4 

8.636 x 10 3 

12 500 

8.705 x 10 16 

2.556 x 10 13 

2.934 x 10 12 


4.436 x 10 4 

1.438 x 10 4 

13 000 

1.091 x 10 17 

3.249 x 10 13 

4.133 x 10 12 


6.249 x 10 4 

2.210 x 10 4 

13 500 

1.315 x 10 17 

3.871 x 10 13 

5.417 X 10 12 


8.191 x 10 4 

3.145 x 10 4 

14 000 

1.525 x 10 17 

4.328 x 10 13 

6.616 x 10 12 

53.2 

1.000 x 10 5 

4.147 x 10 4 


a Emission coefficient for 4300 Arl. 

^Continuum emission coefficient at 4315 angstroms using the formula of 


reference 22. 
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